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Selected problems in heavy quark physics are discussed. The wealth of 
research problems in this field of physics is stressed. 
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Heavy quark physics is a broad and active field of particle physics. Within 
it, hundreds of theoretical papers are produced every year and the produc- 
tion rate keeps increasing. In this short presentation I shall concentrate on 
SpH, recently obtained insights and on open problems. The experimental data 

Qh! quoted without giving the source are either from the 1994 Tables of the 

Particle Data Group [[jj], or from the EPS Conference held in Brussels this 
summer. 

According to the standard model there are six kinds of quarks. In order 
of increasing mass they are denoted u, d, s, c, b, t. The last three are 
considered heavy, because their masses are much larger than Aqcd, i.e. than 
about 0.5 GeV. The mass of the c-quark can be roughly estimated as half 
the mass of the lightest cc quarkonium, which gives m c « 1.5 GeV. This 
in fact it not very heavy — only about three times Aqcd- The mass of 
the recently discovered t quark is m t = (180 ± 12) GeV, which implies that 
the i-quark decays, usually into a VK-boson and a 6-quark, before it has 
time to hadronize. Consequently, the physics of the i-quarks is already well 
understood. The mass spectrum of the heavy quarks causes that most of the 
new ideas apply best to 6-quarks. For t-quarks, the problems are fewer and 



X 



1 Also at the Institute of Nuclear Physics, Cracow. 
This work was partly supported by the KBN grant 2P 302 076 07 



1 



they can be usually solved without making controversial assumptions. For c 
quarks, we are too far from the heavy quark limit, where all the quantities of 
order Aq C d can be neglected compared to the mass of the heavy quark. On 
the long run this may make the physics of the c-quarks more interesting than 
the physics of the 6-quarks, but for the moment it is often just too difficult. 

Let us begin by considering the problem: what is meant by the quark 
mass? 



2 Quark masses 

The standard definition of mass, m = \/E' 2 — p 2 is not applicable to quarks, 
because the energy E and the momentum p on its left-hand-side should be 
measured for free particles. Looking for a free quark not interacting with 
other couloured objects is like looking for one end of a string not attached to 
another end. There is no chance for success. The next choice, when trying 
to define the quark mass, is to make use of the mass parameter mo from the 
Lagrangian. This, however, has corrections. The fact that the quark is part 
of the time a quark-gluon system (the contribution of the gluon loop on the 
quark line) changes the mass by 



s(1) = ma a {n) 



n 



- 7 + log 4vr + log ^- + - 

m 2 3 



(1) 



where 7 = 0.5772 ... is Euler's constant. From this formula one sees two 
difficulties; moreover, there are two others not directly visible. 

• The limit e — > +0 should be taken, thus the formula as it stands does 
not make sense. 

• The scale parameter n is arbitrary. 

• The formula has been obtained using dimensional regularization. There 
are many other methods of regularizing (various cut off procedures, 
putting the theory on a lattice etc.), which yield different formulae. 

• This correction is only the first term of an infinite series, in general 
convergence problems are expected. 



2 



The infinity is eliminated by replacing the mass mo+E' 1 ) by the obviously 
equal number (mo + 5 m) + (E^ — 5m). The trick is to choose 5m so that 
it cancels the infinity in E^. Since nothing is known about m , one can 
assume that E^ does not introduce an infinity in the first term. This recipe 
leaves much freedom in the choice of 5m. Choosing 5m = - one gets the so 
called minimal subtraction mass. Including in 5m also —7 + log(47r), which 
is convenient, one obtains the very popular m mass known as the MS-bar 
mass. Choosing 5m = E^ 1 ) one obtains the pole mass m p etc. Each of these 
masses depends on the scale fi. This arbitrary scale is usually chosen of the 
order of the mass of the quark being considered. For instance, the Particle 
Data Group |]TJ tabulates the quark masses m(m). The differences between 
the various masses are significant. For instance, using the formula for E^ 
one finds for quark Q 

m Q {m Q )=m Q \l — — I. (2) 

Typical values of a s (rriQ) for the heavy quarks are 0.35, 0.20, 0.10 for the 
c, b, t quarks respectively. This gives in the present (very crude) approxi- 
mation the differences between the pole masses and the MS-bar masses 0.17 
GeV, 0.34 GeV and 7 GeV. More careful calculations give for the c and b 
quarks 0.26 GeV and 0.51 GeV ||, while typical values for the i-quark are 
(8 — 9) GeV. An obvious question is: what is the mass found in Fermilab 
for the t quark? The description of the measurement provides an unam- 
biguous operational definition of this mass, but to which of the theoretical 
mass parameters does it correspond? Somewhat surprisingly this problem is 
still controversial. The pole mass, however, seems to be the most popular 
interpretation. 

For the other normalization schemes it is possible to perform analogous 
analyses, therefore, the existence of various renormalization schemes is not a 
serious difficulty 

The convergence problem, however, has been recently found to introduce 
an interesting complication. References can be traced starting from the recent 
review ||. One finds (if one uses dimensional regularization) that the series 
used to define the pole mass is divergent. It can be used as an asymptotic 
series, but then it defines the pole mass only approximately, with an error 
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of about 50 MeV. This is the reason why the MS-bar masses are now the 
popular ones for the c and b quarks. For the t quark the situation is different. 
With present experimental uncertainties an additional uncertainty of 50 MeV 
is irrelevant. On the other hand, the relation between the pole mass and the 
MS-bar mass has a much greater uncertainty. The calculations necessary to 
reduce this uncertainty are possible, but so hard that they have not yet been 
done and are unlikely to be performed it the nearest future. Therefore, if the 
measured mass is the pole mass, expressing it in terms of the MS-bar mass 
would be an unnecessary loss of precision. 



3 Heavy particles 

By heavy particles we mean here particles containing one or more heavy 
quarks or antiquarks. The best studied case is the nonrelativistic approxi- 
mation for the quarkonia QQ. In particular for bottomonia, it is possible to 
get a very good fit to the masses (averages only for the P-states) below the 
threshold for strong decays, for the leptonic widths and for the dipole transi- 
tion matrix elements. One can use the nonrelativistic Schrodinger equation 
with the simple spherically symmetrical potential 

V (r) = asfr + - + c, (3) 
r 

where a, b, c are constants Q. How to make a relativistic theory is still 
controversial. 

For heavy particles containing light quarks the situation is more difficult, 
because for them the nonrelativistic theory does not make much sense. A 
break through has been the idea to use expansions in the inverse of the heavy 
quark mass. For instance, for the mass of a particle with one heavy quark Q 
one finds 



M H = m Q + A + + H T [Darwin + Spin-orbit + Iterll] . (4) 

2m Q 2m Q m z Q 

The leading term is just the mass of the heavy quark. The term of order 
rriQ, denoted A, is the energy of the light component in the colour-field of 
the heavy quark. The heavy quark is here considered as a static source of 



4 



potential. Note the generality of this formulation. The light component may 
be an antiquark, as in valence models of Qq mesons, a pair of quarks, as 
in the valence models of Qqq barions, or a more complicated combination 
of light quarks, light antiquarks and gluons, as in some more sophisticated 
models. The corrections of order 0{m^) correspond to the kinetic energy of 
the heavy quark and to the Pauli interaction of the magnetic moment of the 
heavy quark with the chromomagnetic field created by the light component. 
The magnetic term is responsible for the hyperfine mass splittings in the mass 
spectra. For instance the difference between the B meson and the B* meson 
is that in the first the spins of the heavy meson and of the light component 
give the resultant spin of the particle equal zero, while in the second this spin 
equals one. One finds 

(a-B) = j {Ml* - M 2 B ) « 0.37GW 2 .. (5) 

Since this average should not depend on the mass of the heavy quark, one 
expects a similar value for the (D, D*) system. In fact the experimental 
number is 0.41 GeV 2 . This can be formulated differently: the experimental 
fact that the hyperfine splitting for Q = b is about three times smaller than 
the hyperfine splitting for Q = c, is explained here as a consequence of the 
fact that the c-quark is about three times lighter than the 6-quark. The 
kinetic energy term has no such direct connection to experimental data and, 
therefore, its value is controversial. It can be shown that (p 2 ) > (a ■ B) (|| 
and references given there) and typical estimates are between this lower limit 
and its double. For the terms of order tjiq 2 we have given only the names. 
The first two, the Darwin term and the spin-orbit interaction, are familiar 
from the Dirac theory of the hydrogen atom. The third term is the second 
perturbative iteration of the 0{m~Q L ) term. 

One can apply this approach also to higher resonances. When the light 
component consists of a light antiquark in a P state, its angular momentum 
can be | or | . The parity is plus. Combining that with the spin | and positive 
parity of the heavy quark, one finds four excited states with spins and parities: 
+ , 1 + , 1 + , 2 + . Experimentally one finds two charmed mesons D]* and D%* 
with masses (2423 ±3) MeV and (2458 ±2) MeV respectively and one bottom 
meson B** with mass (5733 ± 17) MeV. A D** meson decays into a pion and 
a D or D* meson. Using angular momentum and parity conservation, as well 
as the information that the pion is produced from the light component, one 
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can see that the mesons with the angular momentum of the light component 
equal |, decay producing a pion in an S-state. Such mesons are broad and 
difficult to observe. The D** mesons with the angular momentum of the light 
component equal §, on the other hand, produce pions in D-states and are 
narrow, because of the suppression of the decay probability by the angular 
momentum barrier. This explains, why only two D** mesons have been 
observed. The hyperfine splitting between these mesons is about 30 MeV. 
Since this is an effect of order ©(m^ 1 ), the corresponding splitting for the B** 
mesons is expected to be about 10 MeV, and indeed cannot be seen at the 
present resolution of 17 MeV. This explains, why for the moment only one 
B** meson has been seen. One also can predict that in order to distinguish 
the two B** mesons, the resolution will have to be improved by about a factor 
of two. 

4 Decays of heavy particles 

Decays of heavy particles are an important source of information about the 
elements of the Cabibbo-Kobayashi-Masakawa (CKM) matrix. From the 
point of view of the standard model these matrix elements are coupling con- 
stants (not all independent from each other!) as fundamental as e.g. the 
electron charge. Where these constants are known, comparison of the the- 
oretical predictions with experiment yields interesting tests of the standard 
model. 

Let us consider the semileptonic decay B — > D*eV. In this decay the 
6-quark, with a probability amplitude proportional to the CKM matrix el- 
ement V c b, goes over into a c-quark. In the process it emits a virtual W~ 
intermediate boson, which decays into the e~, V pair. The problem is to 
extract the modulus \V c b\ from the experimental data. 

In the heavy quark limit the heavy mesons B and D* are similar to hydro- 
gen atoms. In each case the heavy quark sits in the middle, like the proton in 
hydrogen, and the light component surrounds it, like the electron cloud sur- 
rounds the proton in the hydrogen atom. The energy and momentum of the 
W-boson are very large on the scale of the momenta of the light components. 
An analogy would be a 1 MeV photon hitting the proton in hydrogen. In 
this situation the heavy "nucleus" behaves as if it were free. It gets ejected 
with large momentum (on the scale of the light stuff) from its original posi- 
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tion. The 6-quark absorbing (or equivalently emitting) the VF-boson changes 
into a c-quark. Note that since the c-quark is very heavy, large momentum 
does not necessarily mean large velocity. This process, however, is not yet 
the process B — * D* . In order to get the probability amplitude for this de- 
cay it is necessary to multiply the probability amplitude for the ejection of 
the heavy quark by the probability amplitude that the light component of 
the original 5-meson will reorganize itself into the light component of the 
recoiling D*-meson. This is given by the overlap of the two corresponding 
wave functions. Thus, omitting the less interesting (known) terms, the decay 
amplitude is 



Here v and v' denote the initial and the final velocities of the heavy quark. 
In the heavy quark limit these velocities are equal to the velocities of the 
corresponding mesons. The argument uj = v^v'^, which can be interpreted as 
the Lorentz factor of the D* as seen in the rest frame of the B, is a measure 
of the recoil velocity. The overlap factor, known as the Isgur-Wise function, 
is 



Note that the overlapping wave functions of the light components differ only 
by the velocities of their centres. The change of the 6-quark into a c-quark and 
the change of the relative spin orientation of the heavy and light quarks from 
antiparallel to parallel have in leading order no effect on the wave function of 
the light component. The remaining difficulty is how to extract from the data 
the factor \V C b\ without using a specific model for the Isgur-Wise function. 

Two solutions to this problem have been proposed. In the exclusive ap- 
proach one notices that for v — v 1 the two overlapping functions are identical 
and that consequently F(l) = 1 from the normalization of the wave func- 
tion. In this approach one obtains from the data the product | V^F^) and 
extrapolates it to zero recoil, where F(u = 1) = 1. In the inclusive approach, 
one gives up the constraint that the final charmed state must be a D* meson. 
Then the Isgur-Wise function is replaced by the probability amplitude that 
the light component will reorganize itself into anything, which is, of course, 
equal one. Thus, one uses data for the inclusive process B — > X c eV. Here 




(6) 




(7) 
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X c denotes any state containing the quark c. Since the 6-quarks decay al- 
most always into c-quarks, X c can in practice be replaced by X meaning 
anything. We have presented here only the leading term analysis. In prac- 
tice one includes various corrections, which are still somewhat controversial. 
Fortunately they change the calculated values of \V c t>\ by only a few percent. 
Incidentally, the analogous problem of extracting the CKM matrix element 
\Vub\ from the data is much harder and is an active subject of research. 

Let us mention two open problems connected with inclusive decays (cf. 
e.g. ||). Theoretically one finds that the life times of the heavy particles 
containing single 6-quarks are well estimated using the spectator model, i.e. 
neglecting the effect of the light components on the life times. This corre- 
sponds to equal life times for all such particles. It is possible to calculate 
corrections to this result and they turn out to be of a few percent. This agrees 
well with experiment for meson decays, but for the experimental life time 
is only (0.72 ±0.06) of the 6-quark life time inferred from meson decays. The 
theoretical expectation for this ratio is below one, but almost surely above 
0.9. The second problem is the measured fraction of the B mesons, which 
decay semileptonicaly. Theory can reproduce it, but at the condition that 
a large fraction of these decays leads to cc pairs. The average number of c 
and c quarks per decay is experimentally (1.13 ± 0.05), while the theoretical 
number necessary to get agreement with the semileptonic branching ratio 
is 1.3. This difference may seem small, but it should be kept in mind that 
one c-quark is present in almost every 6-decay. Thus what counts is the sur- 
plus over this number. Here the experimental number is less than half the 
theoretical one. 

Finally let us mention the so called rare decays, i.e. the decays, where 
the 6-quark goes over into an s-quark and a photon, or lepton pair. Here 
the theory involves pingwin diagrams, is quite complicated and is still being 
refined, but what is important is that it agrees well with experiment. This 
eliminates many ideas concerning "new physics" i.e. physics beyond the 
standard model. 

5 Production of heavy particles 

Heavy particle production is a broad and active subject. Here we shall only 
mention a few problems, which now are attracting particular interest. 
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The calculated cross-section for the process pp — > ttX at the Tevatron 
is somewhat lower than measured. Since the experimental uncertainties are 
large, however, and since the discrepancy decreases as data improve, this 
does not seem to be a serious problem. 

The ratio of the decay probability of Z° into bb to the decay probability 
of Z° into any hadrons should be about 0.2, because there are five kinds 
of quarks into which a Z° can decay and they all have masses negligible 
compared to the Z° mass. Experimentally 



in agreement with this crude estimate. Precise calculations, however, give 
Rb = 0.2155, i.e. a ratio smaller by about three standard deviations than 
the experimental one. This is considered as a possible problem for the stan- 
dard model. It is interesting that supersymmetry can increase the predicted 
Rb so that it becomes lower than the experimental value by only about one 
standard deviation. If this is the correct explanation of the discrepancy, the 
lightest supersymmetric particles should have masses below 100 GeV and 
there is a good chance of discovering them in the upgraded LEP acceler- 
ator. This is, of course, a bold speculation, but it has recently triggered 
much discussion. Incidentally, the corresponding ratio R c = 0.154 ± 0.07, 
to be compared with the theoretical prediction 0.172. Here, however, the 
experiment is very difficult and a modification of the theory is not plausible, 
therefore this discrepancy is expected to disappear, when data improves. 

Finally let us mention the production of charmonia at the Tevatron. Ac- 
cording to theory those charmonia, which are not decay products of particles 
containing 6-quarks, should be mostly produced in gluon-gluon interactions. 
Such interactions are much more likely to produce P-wave charmonia te- 
states) than S-wave charmonia (^-states). Therefore, the prediction was 
that the direct production of ^-states will be small and that a large majority 
of such charmonia will come from decays of ^-states. Experimentally it seems 
that the direct production of ^-charmonia is much stronger than expected, 
sometimes stronger by more than an order of magnitude. One way out of 
this difficulty is to assume that the cc systems in octet colour states are an 
important intermediate state. 



Rb 




0.2205 ±0.0016 



(8) 



T(Z° — > hadrons) 
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